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We analyze the evolution of the superconducting gap structure in strongly hole doped 
Bai_i,Ki,Fe2As2 between x = 1 and x ~ 0.4 (optimal doping). In the latter case, the pairing 
state is most likely s±, with different gap signs on hole and electron pockets, but with the same 
signs of the gap on the two F-centered hole pockets (a ++ state on hole pockets). In a pure KFe2As2 
(a: = 1), which has only hole pockets, laser ARPES data suggested another s± state, in which the 
gap changes sign between hole pockets (a H — state). We analyze how H — h gap transforms into a 
H — gap as s — >■ 1. We found that this transformation occurs via an intermediate s + is, state in 
which the gaps on the two hole pockets differ in phase by </>, which gradually involves from <j} = ir 
(the H — state) to </> = (the H — h state). This state breaks time-reversal symmetry and has huge 
potential for applications. We compute the dispersion of collective excitations and show that two 
different Leggett-type phase modes soften at the two end points of TRSB state. 



I. INTRODUCTION 

The high interest in iron based superconductors 
(FeSC) is primarily due to two key reasons. The first 
is a hope that the analysis of FeSCs will not only resolve 
the pairing mechanism in these systems but also provide 
important insights into the electronic pairing in a generic 
high-Tc superconductor. The second is a hope to explore 
multi-band structure of FeSCs and discover novel exotic 
superconducting states which have not been observed in 
other systems. Out of such novel superconducting states, 
the most searched for are the ones which break time- 
reversal symmetry. A spin-triplet time-reversal symme- 
try broken (TRSB) ± ipy state has likely been found 
in Sr2RuO^; the spin-singlet d + id TRSB state has 
not yet been observed experimentally, although it was 
once proposed as a candidate state for high Tc cuprate 
superconductors^, and was recently predicted theoreti- 
cally to occur for fermions on hexagonal and honeycomb 
lattices near van-Hove doping^.. 

Several groups already searched for TRSB state 
in FeSCs by exploring the idea that at least in 
some FeSCs both s— wave and d—wave channels are 
attractive^"— liSrii, and that one can, in principle, trans- 
form from s— wave to d— wave pairing by varying system 
parameters - electror^i^ or hole^ doping, hybridization 
between electron pocketsi^, or degree of magnetic scat- 
tering^"*. In between, there is a co-existence regime in 
which both s and d order parameters are present, with 
relative phase ±|, i.e., the system develops a TRSB s±id 
superconductivity. The majority of proposals for s + id 
state are for electron-doped FeSCs, but up to now a d- 
wave superconductivity has not been found in strongly 
electron-doped Ba(Fei_a;Co2:)2As2 nor in KFe2Se2-type 
systems which contain only electron pockets. 

In this communication, we discuss another possible re- 
alization of TRSB state in FeSCs - a purely s— wave state 
with phase difference between superconducting order 
parameters on different Fermi pockets, which is not a 
multiple of tt. The free energy of such a state is sym- 
metric with respect to — > —cf). This Z2 symmetry 



(which corresponds to time reversal since (f> ^ —cj) im- 
plies A — > A*) is broken when the system spontaneously 
chooses (p or —ip. We label such a state as s -I- is. The 
s + is state has been discussed in Refs [l5l - |23l as a generic 
possibility of the superconducting order in the case when 
there are more than two Fermi pockets and as a surface 
state in a two-band superconductor—. We show below 
that TRSB s + is state with varying (f) can be realized in 
strongly hole-doped Bai_j;Ka;Fe2As2 near x = 1. 

We begin by listing several facts about 
Bai_a;K2;Fe2As2. (i) Near optimal doping, x ~ 0.4, 
ARPE S^^'^^ , neutron scattering^!, penetration depth2^ 
and thermal conductivit y^^ ' measurements give strong 
evidence for nodeless, near-constant s± gap, which 
changes sign between hole and electron pockets. This 
is consistent with theoretical calculations^"— 1 ^ ^ 1 , (ii) 
Recent measurements on Bai_2;K2.Fe2As2 with x = 1 
(Refs. and x = 0.93 and x = 0.88 (Ref. M) 

indicate that superconducting Tc most likely remains 
non-zero from x = 0.4 — > 1. (iii) For the x = 1 material 
KFe2As2, ARPES measurements^i^ show that only 
hole pockets are present. According to theory, in this 
situation, both d— wave and s— wave pairing amplitudes 
are attractive^"— and which state wins depends on 
delicate interplay between system parameters, d— wave 
gap is the largest on the hole pocket, which in the 
unfolded Brillouin zone is centered at (tt, tt) (Refs flollTll ). 
and s— wave gap is the largest on the two P— centered 
hole pockets (GCP's), and changes sign between them^. 
The existing experiments point to either c?— wave and 
s— wave gap symmetry: thermal conductivity^!^ and 
specific heaii^ data on KFe2As2 have been interpreted 
in favor of d— wave gap symmetry, while laser ARPES 
measurements^ and other thermal conductivity datcfc^ 
have been interpreted as evidence for s— wave. 

If the gap in KFe2As2 is d— wave, one should obvi- 
ously expect a transition from d— wave to s± state in 
Bai_a;K;j,Fe2As2 as x decreases from 1, and the region of 
an intermediate s + id state at low T^. In this work we 
consider what happens if the gap in KFe2As2 is s— wave. 
At a first glance, one might expect a gradual evolution 
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of the gap structure with x as the symmetry at a; = 1 
is the same as at optimal doping. On a more careful 
look, however, we note that at optimal doping the gaps 
on the two GCP's have equal signs (a ++ state), while 
in s— wave state of KFe2As2 they are of opposite signs (a 
H — state). The issue then is how a H — gap transforms 
into a ++ gap between x = \ and optimal doping. We 
show that this transformation occurs via an intermediate 
s + is state in which the relative phase 4> of the super- 
conducting order parameters on the two GCP's gradually 
evolves between tt (the H — state) and (the +-1- state). 
The system spontaneously chooses either clock-wise or 
counter-clockwise evolution (i.e., positive or negative 0) 
and by this breaks time-reversal symmetry. 

To illustrate the emergence of the s + is state wc first 
consider in Sec. 2 the minimal model with two identical 
GCP's and two electron pockets, all with the same den- 
sity of states A'o, and with the two angle- independent 
repulsive interactions - Uhh between the two GCP's and 
Uhe between hole and electron pockets. A three-band ver- 
sion of this model has been has been considered in Refs. 
Ellll-illlll). The interaction Uhh gives rise to -I — gaps 
on the two GCP's, while Uhe gives rise to an s± state 
with different signs of the gaps on the two hole pock- 
ets. We model the doping dependence by varying the 
strength of hole-electron coupling Uhe and analyze the 
system evolution with Uhe/Uhh- We show that it occurs 
via a TRSB state. In Sec. 3 we extend the model and in- 
clude intra-pocket repulsions and anisotropy between the 
two hole pockets. We show that the TRSB state still ex- 
ists in a certain parameter range, but for non-equivalent 
hole pockets the region of TRSB state is separated from 
Te line. We present our conclusions in Sec. 4. Technical 
details of our analysis are presented in Appendixes A-C. 
In Appendix C we also discuss plasmon mode in a clean 
3D superconductor. 



II. TRSB IN THE MINIMAL MODEL 



Ag is obtained straightforwardly and reads 
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where Uhe = UheNo,Uhh = UhhNo, Nq is the density of 
states, L = In ^nd A is the upper cut-off for the 

pairing. This set can be easily solved. For Uhe > Uhh/V^, 
the eigenfunction with the largest eigenvalue is the H — h 

solution (1,1,-7), where j= ^ + ^ 1 + (j^)" and 

for Uhe < Uhh/V^, is a H — solution (1,-1,0). Precisely 
at Uhe = Uhh/V^ the two states become degenerate and 
0(1,1,-7) + b{l,—l,0) with arbitrary ratio of a/6 be- 
comes an eigenfunction. To see what happens immedi- 
ately below Te at this critical Uhe/uhh we expand the 
Free energy in powers of A^. and Ag. to fourth order 
and obtain (see Appendix A) 



T = To-Ko (lap -I- l&p) + (|a|2 + \b\^y 



~K2\a^ + b^\^ +K3\a\^ 



(2) 



where Ko oc Te - T, Ki^2 > 0, and > K3 > -2K2- 
Minimizing with respect to a and h we immediately ob- 
tain h = ztia\ 
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2K2 



the 



and 



states co- 
exist with relative phase ±7r/2. As a consequence, im- 
mediately below the degeneracy point, the system selects 
an s + is state, which breaks time reversal symmetry (a 
TRSB state). 

Inside the TRSB state we can set Ag to be real and 
Ah, = Ae*'^/2,A^, = Ae-^-^/^^ We solved the set of 
three non-linear gap equations at T = (see Appendix 
B) and found that TRSB state exists between u™" = 
and uj^g'*^ « (1 + HMl iog2). At the lower boundary, 
the TRSB state borders H — state and the relative phase 
reaches = tt, at the upper boundary the TRSB state 
borders state and 6^0. In between. 



The Hamiltonian of the 

— H = Hkin + Hint, 

Uhhhl^^be^k + E»j Uhebl^k^f^k + h 



2 X/fc,Q,/3 



minimal 
where Hkm 
and Hint 



model 



where 
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bxk = J2k^ktaXkii3 and x € {ci, C2, /i, /2}; and 
i,j — 1,2 number the hole pockets (c) and electron 
pockets (/). We define superconducting gaps on two hole 
pockets as Ah-, and A/j^ and the gap on electron pockets 
as Ae, and Ag^. We neglect the angular dependence of 



Uhe in which case Ae 



Ae2 because Uhe for pockets 



ei and 62 are equivalent due to C4 symmetry of the 
underlying lattice. The equivalence between Ae^ and 
Ae2 persists even if we include intra-pocket interactions 
and inter-pocket interaction between the two electron 
pockets. 

The set of linearized equations for Ah, , Ah2 , and 



= ±2 arccos 



Uhh 



,(2«L-«hh)/(2"L"'ih) 



2u 



he 
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We show the evolution of the relative phase 4> on the two 
hole pockets with Uhe/uhh in Fig. [5] 

Combining the results at Te and at T = 0, we obtain 
the phase diagram shown in Fig. [l](a). The TRSB state 
exists in the 'triangle' which begins as a point at £^nd 
extends to a finite interval at T = 0. 



A. Collective modes 

The existence of phase transitions at the boundaries of 
the TRSB state implies that there must be soft collective 
excitations. In a generic multi-gap superconductor there 
are three types of collective excitations: (i) variation of 
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FIG. 1: Qualitative phase diagram for Bai_iKa;Fe2As2 at a; < 1. We model the doping dependence by varying the ratio 
of inter-pocket electron-hole and hole- hole interactions Uhe/uhh which roughly scales as 1 — x. The H — state has gaps of 
opposite signs on the two GCP's and no gap on electron pockets, the -I— I- state is an ordinary s± state in which the gaps have 
opposite signs on hole and electron pockets, and between them is the TRSB state. The gap structures are pictorially presented 
inside each region by vectors placed inside the circles. The magnitudes of the vectors represent |Ai| and the angles represent 
the phases. Cases (a) and (6) are for equal and non-equal intra-pocket interactions {uh^ and Uhj) for the two hole pockets, 
respectively . For (a), the TRSB state starts right at Tc and extends into a finite range at T = 0. For (&), the TRSB region 
splits off from the Tc line and is only accessible at lower temperatures, while immediately below Tc the H — state gradually 
evolves into the -I— I- state as Uhe/uhh increases. 
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FIG. 2; Variation of the relative phase <j) of the gaps on two 
hole pockets with Uhe- This phase is zero for Uhe > w™^, but 
becomes non-zero at smaller Uhe and eventually reaches (p — 
±7r at Uhe = 0. When \(j>\ is between and n, it can be either 
positive or negative, and the choice breaks Z2 time-reversal 
symmetry. The width of the TRSB region is controlled by 
inter-pocket hole-hole interaction u^h and increases when u^h 
gets larger. 

the overall phase, (ii) variations of relative phases of dif- 
ferent gaps (Leggett modes^), and (iii) variations of the 
gap magnitudes. The overall phase mode is coupled by 
long-range Coulomb repulsion to density variations and 
becomes a plasmon^i^. The other modes do not couple 
to density variations are generally either overdamped or 
have energy close to 2A. However, near the boundaries 
of the TRSB state, some of these modes soften. 

We analyzed the dispersion of collective excitations 




FIG. 3: The diagrammatic representation of the equations for 
dispersion of collective modes. The equations for other SAj 
are similar to the one for dAh^ and are not shown. Wavy 
lines - interactions Uij, chain-saw line - Coulomb interaction 
Vq. The bare vertices are not shown. 



in our model by introducing small perturbations in the 
form of pairing and density vertices virith non-zero exter- 
nal momentum and frequency {SA^^ , SA^.-^ , SA^, and Spi, 
i = 1,2,3) and calculating the fully renormalized ver- 
tices (see Fig. [3]). Each SA is generally a complex func- 
tion SAi — Sf + iSf , so for arbitrary momentum q, the 
problem reduces to solving the set of nine coupled equa- 
tions for (5f , (5f, and Spi. We verified, however, that at 
small q, when short-range interactions Uhh, u^e can be ne- 
glected compared to the static Coulomb interaction V{q), 
all three 5pi are equivalent, because the Coulomb repul- 
sion does not distinguish between the different fermions 
(Refs. I4lll43 ). In this approximation, i.e., 5pi = 6p, and 
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the number of equations reduces to seven. 

The equation for <5A/ij^ is graphicahy shown in Fig. [31 



J 



Other equations are similar. In exphcit form we have 



2SP 



3 

3 



(4) 



where 5(0) are bare pairing and density vertices which we 
introduced as small corrections to the Hamiltonian (see 
Appendix C), V{q) is long-range Coulomb potential, and 
the components of the matrix Uij are 



Uhh 2u/ie 
Uhh 2uhe 
Uhe Uhe 



(5) 



Further, H^^ = ^t{q,^) 

^ / (Pkduj /{2tt fTr [gdk,UJ)(J''g^{k + q,UJ + ^)(J^], 
where ct" are Pauli matrices, Qi are Nambu Green's 
function of a superconductor, i S {ci, C2, e}. 

In explicit form we have (see Appendix C for details) 
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A, 2 
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3 
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:Xt 



(6) 



where = In and Xf = 

It is intuitive to reexpress Eq. 2] as 



2A, 



2 E("-').,'5| - 2 E("-').,5;(o) + E n^A' 

3 3 b 

25p = -J: 2NoViq) [n-5f - n-5j + u^PpI (7) 
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FIG. 4: Doping evolution of the frequencies of the relevant 
collective modes at g = 0. The plasmon mode frequency van- 
ishes at g = for all Uhe- Mode Mi describes relative phase 
fluctuation of the two hole gaps. It softens^i2ii2^ at the tran- 
sition point between H — h and TRSB states (at Uhe ~ Uhlit^ , 
indicated by the arrow). Mode M2 describes coupled antisym- 
metric phase fluctuation of the two hole gaps and longitudinal 
fluctuation of the electron gap. This mode lies below twice 
the energy of the electron gap and softens at the boundary 
between TRSB and H — states, at uj^^" = 0. Numerically, the 
energy of the M2 mode becomes small already for Uhe < u™'^ 
because electron gap rapidly decreases with decreasing Uhe- 
The circle represents the case discussed in Ref. [2^ . 



where (5f = [Sf- ,—5l ,5p). This 7x7 set can be cast into 
the form 



Kiq.Vt) S^yT^ (5(0) 



(8) 



where 5 is a 7-component vector with elements 
Sf, —S(, 6p (5(0) is the bare vertex). 
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and 
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Here M = + n33^^ + n33,^ + 2n33. 



coupled to fluctuations of the electron density 5p as 



The dispersions of seven collective excitations are ob- 
tained from the condition T}etK_{q, ft) — 0. 

To adequately describe the full spectrum of all long- 
wavelength collective modes, one has to expand in 
vpq/^, but allow frequency to be of order of A (see Ref. 
|46| and Appendix C). Our goal, however, is limited: we 
want to find the plasmon mode in 2D and the modes 
which soften at the boundaries of the TRSB state. All 
these modes are low-energy modes in the long- wavelength 
limit, and to capture them in our approach, it is sufficient 
to use double expansion in vpq/^ and in il/A. To get 
other modes (or resonances) one needs to search for fre- 
quencies around 2 A. 

In the H — h state, = in equilibrium, and and 
describe phase and magnitude fluctuations, respectively. 
One can easily make sure (see Appendix C) that these 
two sets of fluctuations decouple and there are no solu- 
tions for amplitude fluctuations at 51 <C A. 



The three orthogonal phase modes are 6^ 



Si 



Si, 5{ ^S{ + Si + (2/7)<5|, S', =Si + Si- (2/7)<5|, where 

7 = 2uheLo and = "^^+v/"L+iK.. .^ ^he mode Sl 

is gapped everywhere in the -|--|- phase. The mode 5^ 
describes fluctuations of the overall phase. This mode is 



Vpq'' 



1 



21^2 



4A2 



-Si 



+ S]Sp = 




A 



Sp 



(11) 



The corresponding dispersion is a 2D plasmon with fi^j = 

2 2 

{^Ni^Vq + 1). Observe that the plasmon frequency 
remains the same as in the normal stated. In general, 
Vtpi in a superconductor scales with the density of super- 
conducting electrons and is sensitive to disorder—. In 
our case (clean limit), superconducting density coincides 
with the full electron density, hence ^Ipi does not change 
between normal and superconducting states. 

The mode 5^ describes antisymmetric phase fluctua- 
tions of the gaps on the two hole pockets. The conden- 
sation of this mode signals the transition to the TRSB 
state. In the static limit, this mode totally decouples 

we obtained 



,raax\ 
''he ) 



from density fluctuations. Near Uhe — "-/je 
at g = 0, {^Isif = {%V^I^){2/^/uZ^^?{uhe - u'^ 
Not surprisingly, the antisymmetric phase mode soft- 
ens at the transition point into the TRSB state (where 
Uhe = u^he^)- Wc show the behavior of ^.gi in FigH]). To 
properly obtain the dispersion of this mode, one has to do 
more involved calculations as the combinations of S{,Si, 
and (5|, which decouple at a flnite q, are not the same as at 
q = 0. As a result, the dispersions of Leggett-type modes 
generally depend on the Coulomb interactionSMliiii. 
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Inside the TRSB state, phase and amphtude fluctua- 
tions get mixed up, as was noticed in Refs l2lll23l . This is 
easily seen form Eq. [10] as the off-diagonal components 
which connect the real and imaginary parts of the or- 
der parameter fluctuations, are given by 11^^ which are 
proportional to sin ^ (see Eq. ^ and are non-zero once 

The mode which corresponds to the overall phase 
change is now —{Si — 82 ) sin f + ((5i + (^2) ^os | — ^(5|, 

where in the TRSB state 7 = 2(u/ie/u/i^) cos ^, and (p 
is given by Eq. ([3]). This mode decouples from other 
phase and magnitude modes, but again couples to 5p 
and remains a 2D plasmon. We solved for the remain- 
ing modes and found that the mode 5{ — 82, which 
described antisymmetric phase fluctuations of A/i^ and 
A/12) outside the TRSB region and softened at the up- 
per boundary of the TRSB state, acquires a new func- 
tional form inside the TRSB state, and gets gapped, as 
expected. As uue decreases and cj) increases and ap- 
proaches TT, another mode, indicated as the M2 mode 
in FigHl gets soft. This mode is a coupled oscillation 



of 5t} and (5f 



The first describes longitudinal 



fluctuations of the electron gap, which vanishes at the 
lower boundary of TRSB state, the second describes an- 
tisymmetric phase fluctuations of the two hole gaps (for 
(j) = 7r-20 and 4> « 1, A^i -> Ae'^^'^^ « A(i + (^), and 
Ah^ Ae-<i-^"i K A{-i + ^), and 6^ + 6^ = 26^ = 24> 
describes small deviations from the H — state ). The cal- 
culation of this mode requires some extra care because 
electron gap Ag vanishes at the lower boundary of TRSB 
state, and the expansion in r2^/(2Ae)^ is only valid if the 
mode frequency is below 2Ae fRefl45l). Using the for- 
mal expansion in fl, we obtained the frequency of M2 
mode r2M2 — \/3(2Ae), which is outside the applicabil- 
ity limit of the expansion. A more accurate approach 
is to keep n along with Ag, i.e., replace r2^/(2Ae)^ by 
flV(4A2 - n^). This gives = (V3/2)(2Ae), which 
is below the threshold at 2Ae. We also found another 
low-energy mode using the expansion in J7, however its 
energy is above 2Ae even when we keep fl along with Ag. 
This excitation is then inside the continuum and is not a 
true collective mode. 

We emphasize that the vanishing of Ae is a peculiar- 
ity of the minimal model. In a more general model, 
the TRSB state emerges from the modified H — state, 
in which Ae is already non-zero. Then it is completely 
safe to search for soft modes by expanding in fl/Ai. 



III. BEYOND THE MINIMAL MODEL 

We analyzed whether the TRSB state survives in 
more general cases. As a first step, we included intra- 
pocket density-density interactions u^i, Uh2i and u^- Ap- 
plying the same procedure as before, we found that, 
for u/jj = U}^^ , the phase diagram and the behav- 
ior of collective modes remain the same as in Figs. 



[T] and m the only modification is that at T = 
the lower boundary of the TRSB state now shifts to 

The upper boundary be- 



"■he 



comes u 



he 



(A 



''he 



and ul^ 



X 



where x — 



is the point at which 



TRSB state emerges right at Tc. 

When Uh-^ 7^ , the phase diagram changes qualita- 
tively (see Fig. [T]b). Now one of the hole gaps continu- 
ously evolves from negative to positive along the Tc line, 
passing through zero in between (see Appendix A for de- 
tails). The TRSB state still emerges, but at a lower T, 
and survives as long as intra-pocket interactions remain 
small compared to Uhk (see Appendix B). To simplify the 
presentation, we consider the representative case when 
M/i2,Me = and Ufii ^ Uhh to understand the changes 
to the phase diagram. The phase diagram for a generic 
Uhi 7^ is qualitatively the same as in the case we 
considered. 



We found that TRSB state at T 
ists in an interval between u™" 



Mfefe I 1 

^/2 V 



2/(1 + 6-""! /"''O^ ) 



now ex- 
and M^J^;^^ w 



We also considered anisotropic inter-pocket interac- 
tion Uhh with an extra cos 40 term, consistent with lat- 
tice symmetry^^. This gives rise to cos 46* angular vari- 
ations of Aft,j and A/i^ and may lead to accidental gap 
nodes. The solution of the set of the gap equations for 
Uhi 7^ Uh2 and Uhhid) = Uhhi^ + a(cos6'fti -I- cos 6*^2)) 
is quite involved. However, one can show quite gener- 
ally that TRSB state is confined to low temperatures 
and is separated from the Tc line, like we previously had 
for angle-independent interactions. Immediately below 
Tc, the H — state gradually evolves into ++ state, how- 
ever, now only the average value of the "minus" gap goes 
through zero at some intermediate Uhe, while the gap 
itself does not vanish and just oscillates along the corre- 
sponding pocket. We illustrate this in Fig[SJ 

Inside the TRSB state at T < Tc, the number 
of coupled gap equations equals to nine because in 



general Ah. — Ai 



cos 46*1 ) 



= 1,2. 



= A 



'12 



For — Uh2, we find that A/, 

Ae^'l'/^ {I + {rae-''l' + n) cos For (the 

state), accidental nodes exist if |ra -I- rh| > 1, for </) = tt 
(the -I-- state), they exist li \ra - rb\ > I. In the TRSB 
state, however, |A;iJ doesn't crosses zero and can only 
have gap minima. We illustrate this behavior in Fig. [S] 
for the experimentally relevant case when H — state is 
nodal and state has a full gap. Observe that the 
distance between deep minima gets larger upon entering 
the TRSB state. This behavior is consistent with recent 
laser ARPES studies of doped Bai_^K^Fe2As2 (Ref'32'). 
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s-wave solution 
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FIG. 5: The H — to ++ transition at Tc, with increasing 
Uhe for the case when the two hole pockets are not equiva- 
lent and the interaction u^h has cos 4^ angular dependence. 
The solutions of the linearized gap equations are shown for 
Uhe = 0.02,Mhe = 0.04,u/ie = 0.06, Uhe = 0.09 from left to 
right and top to bottom, respectively. Other parameters are 
a = 0.05,Mhi — 0.2,Uh2 = 0.26, Uhh = 0.2. Note how one 
of the hole gap's average value gets smaller as Uhe increases, 
goes through zero, and re-appears with the opposite sign and 
with small angular variation at larger Uhe- We expect such 
behavior immediately below Tc line in Bai_a;Ka;Fe2 As2, as x 
decreases from one. 




ets, but the same sign of the gap on the hole pockets (a 
++ state in our terminology). In pure KFe2As2 (a; = 1), 
which has only hole pockets, there are experimental and 
theoretical arguments for both d— wave and s— wave gap, 
the latter changes sign between the two GCP's (a H — 
state). We assumed s-wave gap symmetry for KFe2As2, 
consistent with the laser ARPES data?^. The issue we 
addressed is how a -|--f gap on the GCP's transforms 
into a H — gap as a; — 1. We found that, for identical 
GCP's, there is critical point along Tc line at which the 
system jumps form H — to -|— I- state (see Fig. [T]i). At a 
lower T, the transformation occurs via an intermediate 
s + is, state in which the gaps on the two GCP's differ 
in phase by (j) which gradually involves from = tt on 
one end (the H — state) to = on the other end (the 
-|--f state). The system spontaneously chooses either 
or —0 and with this choice breaks time-reversal symme- 
try. We computed the dispersion of collective excitations 
and found that two Leggett-type modes soften at the two 
ends of the TRSB state. We found that the TRSB state 
survives even when the two GCP's are non-identical and 
also when the gap on hole pockets is angle-dependent, 
and even when H — and/or -|--|- states have accidental gap 
nodes. In the former case, near Tc the system gradually 
evolves from the H — to ++ state, but the TRSB state 
still emerges at a lower T (see Fig. ^jp). In the second, 
the nodes get lifted once the system enters into a TRSB 
state (but deep minima remain). The s -\- is state is not 
chiral, but e.g., Kerr effect measurements still should be 
able to detect the breaking of time-reversal symmetry. 
These measurements are clearly called for. 

We acknowledge stimulating conversations with L. 
Benfatto, R. Fernandes, I. Eremin, A. Finkelstein, P. 
Hirschfcld, A. Kamenev, M. Khodas, A. Levchenko, Y. 
Matsuda, I. Mazin, R. Prozorov, J. P. Reid, T. Shibauchi, 
V. Stanev, L. Taillefer, R. Thomale, V. Vakaryuk, and 
M. Vavilov. We are particularly thankful to L. Ben- 
fatto for pointing out an error in our original calcula- 
tion of the collective modes near the lower boundary of 
TRSB state. The research has been supported by DOE 
DE-FG02-ER46900. SM also acknowledges support form 
ICAM-DMR-084415. 
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FIG. 6: The gap evolution in the TRSB state for angle- 
dependent interactions and two equivalent GCP's, in a sit- 
uation when the gap in H — state has accidental nodes and 
the gap in ++ state is nodeless. The nodes disappear once 
the system enters the TRSB state, but deep minima (shown 
by arrows) remain in some range of </) ( or equivalently Uhe)- 



IV. CONCLUSIONS 

We considered the evolution of the superconducting 
gap structure in strongly hole doped Bai_a;Ka;Fe2As2. 
Near optimal doping (x ^ 0.4) the pairing symmetry is 
s±, with different gap sign on hole and electron pock- 



Appendix A: The Free Energy 

We follow a standard procedure and introduce bosonic 
fields Ahi,Ah2, and Ag, which describe fluctuations of 
the superconducting order parameters on the two hole 
and one electron pockets. We decouple four-fcrmion in- 
teractions using a Hubbard-Stratonovic(HS) transforma- 
tion, integrate over fermions, obtain Z = J dAie^'^^^'\ 
and analyze .F[Ai] in the saddle-point approximation. 
For a model with intra-pocket and inter-pocket interac- 
tions within hole pockets (uhnUh2 and Uhh terms, re- 
spectively) and the interaction between hole and electron 
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pockets {uhe term), we obtained 
1 



2uhh - - Uh. 

2(a^,a;;^ + a;;^A;, 

2{uhh - UI-, 



[-2(|A,J2 + |A,J2) 



Uhe 

Uhe ' " 

2{uh^Uh^ -ulf^) 2' 



(a,,a: + a;;^a,) 



A,,,. A* 



which in this case would correspond to either the field (^1 
or ^3.) The condensation of (pi, with 02 = (/"a = brings 
the system into a ++ phase (A/^ — Ah2 — — Ae/7), 
while the condensation of 03 with 02 = 0i = brings 
the system into a H — phase (A^^ — — A/i^, Ag = 0). 
At Uhe = Uhh/V2, ^1 and A2 reach zero at the same 
T, and 0i an 03 condense simultaneously (for this Uhe, 
cos 6 = 1/ VS) . The relative magnitude and the relative 
phase between 0i and 03 are decided by minimizing the 
quartic terms in the Free energy. Plugging in A^ in terms 
of 0i into Eq. lAll neglecting 02, and using Uhe = Uhh/V2 
we obtain 



-2i^|A.p + J G^G^Y.\^-\' (Al) AJ-(4)[0.] = i^i(|0ir + |03H' + if2 



where L = J GG ^ In^, the sum over x runs over two 
hole and two electron pockets, and G — {ioj — e)^^ and 
G={iuj + e)-\ 

Let us first consider the case Uhi — Uh.2 = 0. Then one 
can diagonalize the quadratic part of the Free energy by 
introducing 



4i = cos 9 

62 — sin0- 

A,,- 



Ah^ + Aft., 



A,, 



^h2 



' sin Q Ae 
cos O Ap 



(A2) 



where cos9 = l/^/l + C sinO = C/Vl + C^, and C 
. The action in terms of 



y-hh 



the form 



1 



16< 



AT, 



where 



takes 



Ai|0ip + A2|02p + A3|03|', (A3) 



- 4L 



- 4:L 




Ai = 

A2 = 
A3 = 



(A4) 

Since A2 is strongly negative, the HS transformation 
for 02 does not make sense. Because this field does not 
condense on physics grounds, we just set 02 = (see Ref. 
[49' for more discussion on this). The two other A's change 
sign at some, generally different, temperatures, which de- 
pend on Uhe/uhh- When this happens, either 0i or 03 
condense, depending on whether Ai or A3 changes sign 
first upon lowering T, i.e., increasing L. (This procedure 
is formally equivalent to diagonalizing the linearized gap 
equation to identify the state with the leading eigenvalue 



(A5) 



where Ki ^ ^, K2 = f , K3 = and C > 0. The 

Ki term is isotropic, the K3 term depends on the rela- 
tive magnitudes of 0i and 03 fields, and the K2 term 

if2|02 + 02|2 ^ if2(|0l|^+|03|^ + 2|01P|03|2cOS20). 

depends on the relative magnitude and the relative phase 
between 0i and 03: A positive K2 (our case) selects 
9 = ±7r/2, i.e., if one condensate is real, another is 
purely imaginary. Solving for the amplitudes we find 
103 12 = 101 12 (1 + ^4:3/2X2) = 101 1 V3. The state in which 
both 01 and 03 are present, and the relative phase is not 
or TT is our TRSB state. Eq. IA5I is presented in the 
main text with 0i — >■ a and 03 — ?> b. 

Away from the degeneracy point the quadratic part of 
the free energy takes the form 



J", 



(2)[ 



16 .T , , , 

A + ^ 1011 

3 Uhh, 



A|03P 



(A6) 



where A = 4(1/^^^ — L) and x — 1 — \/2uhe/uhh- The 
leading instability to the left of the degeneracy point (at 
X > 0) is into the 03 state, and to the right of it (at 
X < 0) it is into the 0i state. Once one order sets in, 
it acts against the appearance of the other. Still, we 
found that, e.g., at x > 0, 0i still condenses at Xcr = 
-{16x/3uhh){{Ki + K2)/2K2) = -{lQx/3uhh) * (3/2). 
The corresponding temperature Tcr is smaller than with- 
out K terms, but is still finite. Once 0i becomes non- 
zero, a positive K2 again selects a relative phase of ±7r/2 
between 03 and 0i (which corresponds to the = tt 
boundary for the TRSB state). This consideration leads 
to the phase diagram in Fig. la in the main text. 

We extended this analysis to the case when Uhi = 
Uh2 7^ and found the same results as above. How- 
ever, when Uhi ^ the phase diagram changes qual- 
itatively. To show the new physics and at the same 
time avoid lengthy formulas, we set Uh-^ ^ Uhh] Uh2 = 
and consider Uhe near Uhh/^/2, at which and + — 
phases cross at Tc- Specifically, we set Uh^ — 2yuhh, 
/2)(1 -I- 2c?/), and obtained the phase diagram 



'■he 



"■hh 



to first order in y ^ 1. 

At a non-zero y, the quadratic part of the Free energy 
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reads 



To 



l + y(l-4c)/3 



l + y _ 
Uhh 

1 + y(7- 10c)/3 



^101 



i 102 



V3i 



(05(02- V20i)+ c.c) (A7) 



The 02 rnode is again non-critical, and 02 can be sent to 
zero. For the remaining two modes, we have 



^2[0.] = 4fii±^-L)|03p 

V Uhh J 

l + y(l-4c)/3 . 



4?/ 



Uhh 
03 0^ + +C.c) 



(A8) 



V^Uhh 

Diagonahzing this quadratic form by 

01 = V'l COS »7 + V'3 sin 77 03 = —01 sin 77 + "03 cos 77, (A9) 

we obtain tan 277 = \/3/ (1 + 2c) . Taking the positive root 
tan77 = -ij ^(1 + 2c)2 + 3 - (1 + 2c) , we obtain 



1 + f (2(1 - c) - V(l + 2c)2+3 



Uhh 



1 + ^ 
^ ^ 3 



(2(1 - c) + V(l + 2c)2 + 3) 



(AlO) 



We see that the temperatures at which ipi and 03 
modes condense are different and V'l mode condenses 
first for all values of c. The ipi mode condenses at L^^ — 
{l + Si{c))/uhh,yf^eK Si{c) = (8/3)(l-c-x/l + c + c2), 
and the 1/^3 mode condenses at L^^ = (1 + S3{c))/uhh, 
where S'3(c) = (8/3)(l - c + Vl + c + c^). We plot the 
temperatures at which the prefactors for IV^ip and |'03p 
terms vanish in Fig [7] The condensation of ipi field leads 
to a superconducting state in which all three gaps A^^ , 
A/12 J Ae are generally present and are different from 
each other. At large positive c (i.e., at larger Uhe) the 
state immediately below the condensation temperature 
of 01 is close to the ++ state, with A/i^ A/i^ and Ae 
of opposite sign compared to A/j^ and A/12 . At large 
negative c (smaller Uhe) the state immediately below the 
condensation temperature of "01 is close to the H — state, 
with A/jj w — A/i2 and smaller Ag. In between, the con- 
densed state is a mixture of H — h and H — states. In par- 
ticular, for c = 0, Ae = —Ah2/V^ and A/^ = 0, i.e., the 



gap on the hole pocket, for which we kept intra-pocket 
repulsion, vanishes. We analyzed the form of the con- 
densate for various c (i.e., various Uhe/uhh) and found 
a continuous evolution, in the process of which one of 
hole gaps gets smaller, passes through zero, and then re- 
emerges with the opposite sign. Specifically, we found, 
right below Tc for the "01 mode. 



A,a 



i7i2 



V2 



Ahi - A 



h2 



Ahl 



= (1 + 2c) - v/3+(l + 2c)2 (All) 



Without quartic terms, the modes ipi and ips are de- 
coupled and the system undergoes two superconducting 
transitions at L^, and L^^. The mode which condenses 



at 1/^3 is almost H — h state at large negative c, almost 
H — state at large positive c, and a mixed state in be- 
tween. E.g., at c = 0, the "03 condensate has components 
A,,, = -2Ah,,Ae = AhjV2. 

The situation changes when we include quartic terms 
into consideration. We use Eq. (|A5I) as an input, sub- 
stitute 0i_3 in terms of ipi^s via and obtain T4[ipi]. 
Carrying out the calculations, we find that the four- 
fold term contains a linear piece in -03 in the form 
2if3 sin277COs2 7y|'0ip|i/'3| cos6'i3, where ^13 is a relative 
phase between the condensates of "01 and 4'3- This term 
acts as an "external field" for -03 and makes 7/^3 non-zero 
once -01 condenses. Because K3 < 0, the system initially 
selects 6*13 = 0, i.e., 03 field emerges with the same phase 
as 01. This implies that the state immediately below 
breaks a U(l) gauge symmetry (the overall phase 
gets fixed), but time-reversal symmetry remains unbro- 
ken. The situation changes, however, when the temper- 
ature gets lower and ips grows. The full dependence of 
J^4[7/'i] on 013 is in the form 

J^4 [ipi] — 2K3 cos 013 sin 277 

^IMi^sl (|0'ipcos2r;+ |7/;3|2sin2 77) 

+ cos2 0i3|7^i|2|7/>3|2 {4K2 + K3sm^2Tj) (A12) 

Analyzing this form, we immediately find that the pref- 
actor for cos^ ^13 is necessary positive. Minimizing with 
respect to 0i3, we then find that, at some finite 7/^3, the 
equilibrium value of 0i3 shifts from 0i3 = to a finite 
013 = ±6, 5^0. For large and small c, this happens al- 
ready at small tp^ , which are well within the applicability 
of the expansion in powers of ip- Thus, for large positive 
c, the critical l-ipa] = |7/'i|/(\/3|l + 2c|). 

Once the system selects a non-zero 0i3 , it breaks addi- 
tional Z2 symmetry by selecting either positive or nega- 
tive value of the relative phase 0i3 . The Z2 breaking then 
implies that time-reversal symmetry is broken, i.e., once 
013 becomes non-zero, the system enters into a TRSB 
phase. The region of this phase shrinks as Uh^ increases 
but definitely remains finite as long as Ti^i << Uhh, i-C, 
as long as our parameter y is small. 

When both Uh^ and Uh2 are non-zero, the calculations 
become more involved, but the physics remains the same. 



10 



We also analyzed the effect of adding intra-pocket in- 
teraction Me for electron pockets. Like in the case of 
Uhi — Uh.2, a non-zero shifts the lower boundary of 
the TRSB state to a finite Uhe- There is one new effect 
compared to the case Uh-i — Uh^: because now (when 
M/ii 7^ Uh2), remains non-zero to the left of the lower 
boundary of the TRSB state, the mode which describes 
longitudinal fluctuations of Ag, no longer strongly cou- 
ples to antisymmetric phase fluctuations of the two hole 
gaps, and the mode which softens at the lower boundary 
of TRSB state becomes a pure Leggett-type phase mode. 



:-wave solution 



Equal hole pockets 



Uriequal hole pockets 



-leading 
subleading 



-2 -1 



FIG. 7: Temperatures (T) at which the prefactors for the 
quadratic terms in Ginzburg-Landau expansion for the two 
critical fields change sign. The parameter c measures the de- 
viation of hole-electron interaction Uhe from the critical value 
(= Uhh/ \/2) . Left panel - two equivalent hole pockets {y — 0) . 
In this situation, the condensation of one critical field leads 
to H — order, the condensation of the other leads to -I— f or- 
der. The two lines cross at the critical Uhe- In the presence 
of mode-mode coupling, the emergence of one order tends to 
prevent the emergence of the other, and the actual tempera- 
ture, at which the second order emerges, gets smaller (black 
line). We found (see text) that below the black line the two or- 
ders lock into TRSB state. Right panel - non-equivalent hole 
pockets {y — 1/8). The eigenfunctions reduce to pure H — and 
-|— |- only at large |c|, while in the region labeled 'mixed' the 
system gradually transforms form the H — to H — h order with 
one of the hole gaps going through zero in between. The lines 
at which the prefactors for the quadratic terms vanish now 
do not cross. Due to mode-mode coupling, the order which 
appears first now induces second order, i.e., both are present 
immediately below the actual Tc line, with a relative phase of 
or TT, i.e., time-reversal symmetry is not broken at Tc- The 
TRSB state still emerges, but at a lower T (below the black 
line). 



Appendix B: Non-linear gap equations at T = 

The key goal of the analysis is to show that TRSB 
state, which starts as a point along Tc line, extends to a 
finite range of system parameters at T = 

The set of non- linear gap equations in a generic model 
with inter-pocket interactions Uhh, Uhe, and intra-pocket 
interactions Uhi, Uh2, and Ue is shown diagrammatically 
in Fig [HI Each anomalous vertex is a gap A^;, which, in 
general, is a complex variable {x = ft.1,/12, and e), and 
each fermionic bubble is a sum of normal and anomalous 



Green functions 



2 -h £? ' 



(Bl) 



where Ex — -y/e^ -f | Aj^ , is the fermionic dispersion 
near the pocket x, and Qa^p = icr^^. Evaluating the dia- 
grams, we obtain at T = 



A,, = 
A,, = 
A, = 

where Lx = h 



-Ufii^hiLi - Uhh^h2L2 - 2u/ieAeie 

2u/ieAeie 
UeAeLe (B2) 



-Uhh^hiLi - u/i2A/i2i2 




( ^ 



'hh + 




■'4 



FIG. 8: Diagrammatic representation of the set of non-linear 
equations for the gaps A^^ and Ae^ (viewed as anomalous 
self energies). In our case Ae^ — Aej = A^. The equation 
for Afi2 is similar and not shown. The double headed arrows 
correspond to the anomalous Greens functions; The single and 
double solid lines and the single and double dotted lines are 
anomalous Green functions for fermions near the hole pockets 
(/ii,2) and near electron pockets (£1,2), respectively. 



1. The symmetric case 

Consider first the symmetric case u^-^ — ■ Then 
A/ij = A/12 — ^ and Li = L2 = L. Without loss of 
generality, the overall phase can be set such that Ae is 
real. The two hole gaps must then satisfy A/i^ — A^^ , i.e 
in general Ah^ = Ae*'^/^, Ah^ = Ae''"^!'^. The electron 
gap Ae also scales with A, and we write Ag = —7 A, 
in which case Le = i — Iw^. The three variables A, 7, 
and (/) are the solutions of the set of three non-linear gap 
equations (we recall that L = log^). We have, from 



(Eq. |B2 



[1 - {uhh - Uhi) sin(0/2) 
[1 + [uhh + tihi ) L\ cos(0/2) 

[1 + Meie]7 





lUhelEe 

2uhe cos((/)/2)L 



(B3) 
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For the H — state, = tt, and we have 7 = and 
L — l/{uhh — Uhi)- For the ++ state, = 0, i is ap- 
proximately the smallest positive solution of 

1 + (Ue + U,ih + UhJ)L + {Ueiuhh + "/ij " 4li^e) = 

(B4) 

and 7 is the solution of 7(1 + Ue{L — log^) = 2uheL. 

For the TRSB state, (j) is different from and tt, and 
we have 



L = 



1 



Uhh - Uhi 
Uhh 



7 = 2 



UeUhh 
UheL (f) 

cos — 

1 + UeLe 2 



(B5) 



The upper and lower boundaries of the TRSB state 
are obtained by matching the TRSB solution and the 
solutions for the +- 1- and H — states, respectively. This 
gives u^^^ and u™", which we presented in the main 
text. 



2. Non equivalent hole pockets 

For Uh2, A/i, = Aic"'^*/^, and both Ai^2 and 

01,2 are generally different. The analysis now involves 
five variables (two complex A^. an one real Ag, and is 
quite involved. However, less efforts are needed to just 
prove that TRSB state exists because near its upper and 
lower boundaries (f)i and 02 approach zero or differ by tt, 
respectively, and one can expand in the deviations from 
equilibrium 0i's. 

As an example, consider the system near the upper 
boundary of the TRSB state. Here 0i and 02 are both 
small. Expanding in the set of complex equations (jB2[) 
for and Ag to linear order in 0i,2) and separating 
real and imaginary parts, wc obtain, from the imaginary 
parts, 

Ai (1 + UhlLi) 01 + L2 A202 = 
A2 (1 + Uh2L2) 02 + iiAi0i = 
AiLi0i + A2L202 = (B6) 

Combining, e.g., the first two and the last two equations 
and each time setting the determinant to be zero and 
combining with the third equation in (IB6[) . we immedi- 
ately obtain 



r 1 2A 1 

Li = log-r- = : 

Ai Uhh - Uhi 

T 1 2A 1 

L2 = log— = , 

A2 Uhh - Uh2 



(B7) 



The real parts of the same set of Eas. (jB2p can be 
evaluated at 0i =02 =0. The first two equations of 
the set (IB2I) with real A^,. — Ai are identical for Li 2 



(and Ai,2) given by (jB7[) and using them we can ex- 
press AgLe = Ag log in terms of various couplings u. 
Solving for Ag and substituting the result into the last 
equation in (|B2I) we obtain the expression for Uhe = i*5Je"^ 
for the upper boundary of the TRSB state. The result 
for M™g"^ for Uh2 = Ug = and Uhi << Uhh is presented 
in the main text. The result for the lower boundary of 
the TRSB state, u™" is obtained in a similar manner, 
by expanding near 0i,2 = tt. 



3. TRSB state for angle-dependent interaction 

Our primary interest is to study how the TRSB state 
is modified if outside this state the gaps on the two F- 
centered hole pockets have angular dependence and even 
accidental nodes, if this dependence is strong enough. 

To focus on this physics and avoid lengthy formulas, 
we ignore potential anisotropy of intra-pocket interac- 
tions Uft. and Mg and of electron- hole interaction Uhe^ 
and only include the anisotropy of the interaction Uhh 
between the two F-centered hole pockets. By symme- 
tryii, angle-dependence of Uhh comes in the form 

Uhh{k,p) =it/ih(l + 2acos46'fc-F2Q;cos46'p)-H... (B8) 

where dots stand for cos 85, etc terms which we neglect. 
The most general solution for the hole gaps for this form 
of the interaction is 



Ai (e 



A., 

Ah2 = A2(e'' 
Ag = A3 



+ Tie**" cos 461) 
+ r2e*'^='' cos 461) 



(B9) 



where without loss of generality we can set A^ and 
to be positive. As before, we select Ag to be real by 
adjusting the overall phase. 

To obtain the gaps in the TRSB state for arbitrary in- 
teractions M, one has to solve the set of nine coupled 
equations, which can only be done numerically. One 
can, however, still find an analytical solution for the case 
u/ii — Uh2- In this situation, two hole pockets are equiv- 



alent, and one can easily show that A/^ 



We 



verified that the set of non-linear gap equations is satis- 
fied if we use the following ansatz 



A;,, = A^^ = Ae»*/2 (1 . 
Ag = -7A 



[rae 



Tb) COS 46 



(BIO) 



This ansatz contain five unknowns (A, 7, Ta, fb, 0). Sub- 
stituting these forms into the set of non-linear gap equa- 
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tions (Eq. (IB2p with ujih given by (jB8|) . we obtain 
= -2auhh y ( 1 + Tfc cos 46*) 
-2auhhra / i(?cos46' 



r-6 



cos ■ 



(u/i/i^e +UftjLe(l + r6 cos 40) cos- 
ra {uhhA-e + Uhi)Le COS 4:9 COS - 

2uheLe'-f 

{uhhAe -Uhi)Le{l + n cos 46) 



ra {uhhAe ~ UhjLe cos 40 



7 



2uh, 



J LeCOS^{l + {ra + n)cos49) (Bll) 



where Lg = In 
we have r„ = 



2A 

|A(e)| 



andAg = l+2acos40. When a = 0, 
= 0, and the other three equations 



coincide with what we had in the isotropic case. 

We analyze this set both analytically and numerically, 
and found that TRSB state (the one with different from 
zero or tt) still exists, at T = 0, in some range of w^e, even 
if the hole gaps in H — and/or ++ states have acciden- 
tal nodes. However, in the TRSB state, the gap ampli- 
tude has minima but no nodes, simply because |A/iJ = 
\Ah^\ = A'^{{l + {raCOS(l) + rb)cos4eY +rlsm^ (f>cos'^ 49) 
never hits zero when sin is non-zero. We discuss this in 
the main text. 



Appendix C: Collective modes 

In this Appendix we present some details of the deriva- 
tion of the dispersion of collective modes. 



We consider the minimal model with two equal hole 
pockets and two inter-pocket interactions Uhe and Uhh- 
The extension to more general cases is straightforward, 
but the formulas become more cumbersome. 

We include both the pairing interactions {uhe and Uhh) 
and 2D long-range Coulomb interaction Vq = A2/\q\, 
A2 = 27re^. To obtain the dispersion of collective modes, 
we add to the system a small frequency and momentum- 
dependent perturbation (the bare terms) 



pert 



-t-^[ci O C2] 

fe fe 

+ ^[ciO/i]+^[ciO/2] 
k k 

+6piO)J2(c\ci + ... + .412) (CI) 



where SA, = 5A,((j, r2)e*(^^*-i "") and Ap = 
Sp{q,n)e'^^*~'^''^\ compute fully renormalized 6 A 
and 6p, and obtain collective modes as the poles of the 
generalized susceptibility. Alternatively, the collective 
modes can be compute d by e xtending HS approach to 
finite q and D,, see Refs. I44l46ll49l . 

The field Sp{q, il) = 5p is real, while 6Ai{q, ft) is gen- 
erally complex and it is instructive to split it into real 
and imaginary parts: SAj{q,V,) — 6^^ + iSj. If the 
equilibrium gap Aj is real, (5j^ and 5j describe ampli- 
tude (longitudinal) and phase (transverse) fluctuations 
of the gap. If the equilibrium gap is complex, each of 
and 6j describes amplitude and phase fiuctuations. 

In particular, if in equilibrium Ahi = Ae*'^/^,A/i2 = 
Ae~*'^/^,Ae — — 7A, the relation between (5j^, 6j and 
the changes of the amplitudes and the phases of the three 
gaps [| A;,J ^ A + nih, , </./2 ^ (f>/2 + (fi^, ; 
Ah,^A + nih, , ~(f)/2 -4>/2 + (j,h^ ■ 

I Ae| -> -7A me, ^ 0e] is 
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(C2) 



Each of the bare vertices gets renormalized by the pair- 
ing interactions and long-range Coulomb interaction. At 
weak coupling, only ladder-type particle-particle renor- 
malizations and small q particle-hole renormalizations 
are relevant. Collecting the relevant diagrams (see Fig. 
U] in the main text), we obtain the set of coupled equa- 
tions for fully renormalized vertices SAi — 5^ + iSj and 



Sp — Sp as we said in the main text. 

The seven branches of collective excitations are ob- 
tained from the condition that Deti^(q, ft) — 0. Two of 
these branches are fluctuations of the overall phase and 
of the total density, the others are three longitudinal gap 
fiuctuations and two different fluctuations of the relative 
phases of the three gaps. Some of these fluctuations de- 
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couple from the others, but some are coupled. 

The components of n^^''(g, Q) can be represented in the 
Nambu formalism as 

(C3) 

where ui is the fermionic Matsubara frequency, ct* are the 
Pauli matrices, and 



Gi{k,uj) 
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~F,{k,uj) 
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(C4) 



where 
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Gu 
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It 



luj + ek,i 
iuj - ek,i 
A,. 



Ef 



F^ - -- 



A* 



Ef 



(C5) 



Evaluating the integrals, we find that 21 components of 
n are non-zero. 

To properly describe all collective excitations, one 
should keep the frequency to be of order A, as some of 
the modes exist only as resonances at > 2A. Our 
goal, however, is more focused as we are only interested 
in the 2D plasmon mode and in the modes which soften 
at the boundaries of TRSB state. These modes are the 
solutions of Detjf (g, 57) = at small 57, and to get these 
modes one can safely expand in both vpq/ A and in 51/ A 

Evaluating the integrals and converting from Matsub- 
ara to real frequency axis we obtain the expressions for 
and K_[q, 51) at small 51 and q, which we presented 



in Eqs. ^ and pUj) in the main text. 

Solving for Pet J^fq, 51) — 0, we obtain seven branches 
of collective excitations, which we discuss in the main 
text. One can show quite generally that fluctuations of 
the overall phase and of the total density are coupled 
to each other but decoupled from other five branches of 
collective excitations. One of coupled oscillation of the 
overall phase and the total density is a plasmon mode 
(see the main text). Among the other five modes, longi- 
tudinal and transverse fluctuations decouple in ++ and 
H — phases, but couple in the TRSB state. This coupling 
leads to a peculiar structure of low-energy collective ex- 
citations near the boundaries of the TRSB state. We 
present the results in the main text. 



1. Plasmon mode in a 3D superconductor 



oscillation of a phase of a superconductor order parame- 
ter and an electron density) in a 3D superconductor. In 
3D, plasmon frequency tends to a finite value at g 0, 
and the approximation 52 ^ A, which we used in the pre- 
vious subsection, is not applicable, at least in the clean 
limit. 

In the dirty limit, the plasmon frequency is small (it 
can be much smaller than A). A general gradient expan- 
sion analysis in this case shows^ that the plasma fre- 
quency scales with the density of superconducting elec- 
trons (the "superfluid density"). In a clean limit, su- 
perfluid density coincides with the full density, and it is 
reasonable to expect that the plasma frequency remains 
the same as in the normal state. 

That the plasma frequency is not renormalized in the 
clean limit and at T = has been argued by Anderson 
back in 1958 on general grounds (Ref.-^) and has been 
shown explicitly by Ohashi and Takada using an RPA 
formalism, extended to a superconduting stated. We re- 
produce this result in a direct diagrammatic approach, 
similar to the one we used in the main text for the 2D 
case. For briefness we consider the case of a single-band 
s-wave superconductor. The extension to multi-band sys- 
tems is straightforward. 

We follow the same strategy as in the main text - in- 
troduce bare particle-particle and particle-hole vertices, 
which correspond to small variations of a superconduct- 
ing gap and a total density {5/S. — 5^ + iS^ and Sp, re- 
spectively), and express the full vertices in terms of the 
bare ones, using dimensionless u < for the pairing in- 
teraction and Vq = A-ijq^ for Coulomb interaction in 3D, 
with A3 = 47re^ . The diagrams for the vertices are shown 
in Fig. [i 





o < 



FIG. 9: Diagrammatic representation of the coupled equa- 
tions for fluctuations of the total density &p and the SC order 
parameter 5 A (and 5A*) for the one band case. The solid 
and dotted wavy lines represent the pairing interaction u < 
and unscreened Coulomb interaction Vq . The lines with single 
and double arrows represent the normal(G) and anomalous(F) 
Green functions. The coupling is due to GF terms which are 
non-zero when g , f2 7^ 



For completeness, we also present the diagrammatic 
derivation of the dispersion of a plasmon mode (a coupled 



Like in the previous section, we introduce the vector b 
with the components (5^,-5^, and and write the full 
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vertex S in the same was as in 



but now with 



Eq. IC7l now becomes, to the leading order in q: 




The zeros indicate that the magnitude fluctuations 6^ do 
not couple to the phase and density fluctuations {5^ and 
dp terms). The last two fluctuations, however, couple 
to each other. The dispersion of the collective modes 
are again obtained from the condition Det_K(g, il) = 0. 
The mode which corresponds to coupled phase-density 
oscillations is obtained from 



2 . 

- - 
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Expanding only in q, we get 
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where E' = Ve^ -|- A^, A is the upper cutoff, and Q'^ = 
{(vF.q)^) = in 3D (and ^ in 2D). In H^^ we have 
used the BCS gap equation that tells us 
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Also, 
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we immediately find that 
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n 

2A 



(C13) 



and hence 



(CM) 



which is the same result as in the normal state. Substi- 
tuting the expressions for Vq — Aire'^/q'^, — Vpq'^/3, 
Nq = mpp / and using the relation between pp and 
the density of fermions ^^./(BTr^) — n, we obtain 



Anne 



pi 



(C15) 



which is the same plasma frequency as in the normal 
state. This result is well-known starting from the Ander- 
son work'*^. Like we said, our goal was just to demon- 
strate how this result can be re-derived in a direct dia- 
grammatic approach. 

At a finite T < Tc and/or in the presence of impu- 
rity scattering, coupled density and phase fluctuations 
are more complex, and near Tc there exists a weakly 
damped, near-gapless Carlson-Goldman mode^S. The 
evolution of plasma oscillations with increasing T and/or 
impurity scattering are not fully understood as only the 
cases SI = fipj and $1 ^ A <C have been analyzed in 
detail (see, e.g., Ref. 113). The diagrammatic approach 
which we present here offers the way to obtain the results 
for all T and also with and without impurity scattering. 



^ A.P. Mackenzie, Y. Maeno, Rev. Mod. Phys. 75, 657712 
(2003). 

^ R.B. Laughlin, Phys. Rev. Lett. 80, 5188(1998). 

^ R. Nandkishore, L.S. Levitov, and A.V. Chubukov, Nat. 



Phys. 8, 158 (2012). 
* W.-C. Lee, S.-C. Zhang, and C. Wu, Phys. Rev. Lett. 102, 
217002 (2009). 

^ K. Suzuki, H. Usui, and K. Kuroki, Phys. Rev. B 84, 



15 



144514 (2011). 

^ S. Graser, T.A. Maier, P.J. Hirshfeld, D.J. Scalapino, New 
J. Phys. 11, 025016 (2009). 

F. Wang, H Zhai, Y Ran, A. Vishwanath, and D-H Lee, 
Phys. Rev. Lett. 102, 047005 (2009). 
* R. Thomale, C. Piatt, W. Hanke, and B. Andrei Bernevig, 
Phys. Rev. Lett. 106, 187003 (2011). See also C. Piatt, 
C. Honerkamp, and W. Hanke, New J. Phys. 11, 055058 

(2009) . 

^ R. Thomale, C. Piatt, J. Hu, C. Honerkamp, and B. A. 

Bernevig, Phys. Rev. B 80, 180505(R) (2009). 
^° R. Thomale, C Piatt, W. Hanke, J. Hu, and B. Andrei 

Bernevig, Phys. Rev. Lett. 107, 117001 (2011). 

S. Maiti, M. M. Korshunov, T. A. Maier, P. J. Hirschfeld, 

and A. V. Chubukov, Phys. Rev. B 84, 224505 (2011); 

Phys. Rev. Lett. 107, 147002 (2011). 
^2 C. Piatt, R. Thomale, C. Honerkamp, S.-C Zhang, W. 

Hanke, Phys. Rev. B 85, 180502 (2012). 
" M. Khodas, A. V. Chubukov, Phys. Rev. Lett. 108, 247003 

(2012). 

R. M. Fernandes and A. J. Millis larXiv:1208.3412i (2012). 

V. Stanev and Z. Tesanovic, Phys. Rev. B 81, 134522 

(2010) . 

^® D.F. Agtenberg, V. Barzykin, and L.P. Gorkov, Phys. Rev. 

B 60, 14868 (1999). 
" T.K. Ng and N. Nagaosa, Europhys. Lett. 87, 17003 

(2009). 

Y. Tanaka and T. Yanagisawa, Solid. State. Comm. 150, 
1980 (2010); T. Yanagisawa, Y. Tanakal, I. Hase, and K. 
Yamaji, J. Phys. Soc. Jpn. 81 024712 (2012). 
X. Hu and Z. Wang, Phys. Rev. B 85, 064516 (2011). 
2° S.-Z. Lin and X. Hu, Phys. Rev. Lett. 108, 177005 (2012). 
J. Carlstrom, J. Garaud, and E. Babaev, Phys. Rev. B 84, 
134518 (2011). 

^■^ G. Livanas, A. Aperis, P. Kotetes, G. Varelogiannis, 

I arXiv: 1208.2881 (2012). 

V. Stanev, Phys. Rev. B 85, 174520 (2012). 
2" I. Bobkova and A. Bobkov, Phys. Rev. B 84, 134527 

(2011) . 

H. Ding et al., Europhys. Lett. 83, 47001 (2008). 
K. Nakayama et al., Phys. Rev. B 83, 020501 (2011). 
A.D. Christianson, E.A. Goremychkin, R. Osborn, S. 
Rosenkranz, M.D. Lumsden, CD. Malliakas, L.S. Todorov, 
H. Glaus, D.Y. Chung, M.G. Kanatzidis, R.I. Bewley, and 
T. Guidi, Nature 456, 930 (2008). 

R. Khasanov et al., Phys. Rev. Lett. 102, 187005 (2009). 
X. G. Luo, M. A. Tanatar, J.-Ph. Reid, H. Shakeripour, 
N. Doiron-Leyraud, N. Ni, S. L. Budko, P. C. Canfield, 
Huiqian Luo, Zhaosheng Wang, Hai-Hu Wen, R. Prozorov, 
and Louis Taillefer Phys. Rev. B 80, 140503 (2009). 
^° J.-Ph. Reid, A. Juneau-Fecteau, R. T. Gordon, S. Rene 
de Cotret, N. Doiron-Leyraud, X. G. Luo, H. Shakeripour, 
J. Chang, M. A. Tanatar, H. Kim, R. Prozorov, T. Saito, 

H. Fukazawa, Y. Kohori, K. Kihou, C. H. Lee, A. lyo, H. 
Eisaki, B. Shen, H.-H. Wen, and L. Taillefer, Supercond. 
Sci. Technol. 25, 084013 (2012). 

I. I. Mazin, D. J. Singh, M. D. Johannes, and M. H. Du, 
Phys. Rev. Lett. 101, 057003 (2008); K. Kuroki, S. Onari, 
R. Arita, H. Usui, Y. Tanaka, H. Kontani, and H. Aoki, 
Phys. Rev. Lett. 101, 087004 (2008); A. V. Chubukov, D. 
V. Efremov and I Eremin, Phys. Rev. B 78, 134512 (2008); 
V. Cvetkovic and Z. Tesanovic, Phys. Rev. B 80, 024512 



(2009); J. Zhang, R. Sknepnek, R. M. Fernandes, and J. 
Schmahan, Phys. Rev. B 79, 220502(R) (2009); I.I. Mazin 
and J. Schmahan, Physica C., 469, 614 (2009); A. F. Kem- 
per, T. A. Maier, S. Graser, H-P. Cheng, P. J. Hirschfeld 
and D. J. Scalapino, New J. Phys. 12, 073030 (2010); P. J. 
Hirschfeld, M. M. Korshunov, I. I. Mazin Rep. Prog. Phys. 
74, 124508 (2011); A. V. Chubukov, Annu. Rev. Cond. 
Matt. Phys., 3, 57 (2012). 

K. Okazaki and S. Shin, private communication. 
T. Sato, K. Nakayama, Y. Sekiba, P. Richard, Y.-M. Xu, S. 
Souma, T. Takahashi, G. F. Chen, J. L. Luo, N. L. Wang, 
H. Ding, Phys. Rev. Lett. 103, 047002 (2009). 
K. Okazaki, Y. Ota, Y. Kotani, W. Malaeb, Y. Ishida, T. 
Shimojima,T. Kiss, S. Watanabe, C.-T. Chen, K. Kihou, 
C. H. Lee, A. lyo, H. Eisaki, T. Saito, H. Fukazawa, Y. Ko- 
hori, K. Hashimoto, T. Shibauchi, Y. Matsuda, H. Ikeda, 
H. Miyahara, R. Arita, A. Chainani, S. Shin, Science, 337, 
1314(2012). 

S. Maiti, M. M. Korshunov, and A. V. Chubukov, Phys. 
Rev. B 85, 014511 (2012). 

J.-Ph. Reid, M. A. Tanatar, A. Juneau-Fecteau, R. T. Gor- 
don, S. Rene de Cotret, N. Doiron-Leyraud, T. Saito, H. 
Fukazawa, Y. Kohori, K. Kihou, C. H. Lee, A. lyo, H. 
Eisaki, R. Prozorov, Louis Taillefer, Phys. Rev. Lett. 109, 
087001 (2012) 

A. F. Wang, S. Y. Zhou, X. G. Luo, X. C. Hong, Y. J. 
Yan, J. J. Ying, P. Cheng, G. J. Ye, Z. J. Xiang, S. Y. Li, 
and X. H. Chen, arXiv 1206.2030. 

M. Abdel-Hafiez, V. Grinenko, S. Aswartham, I. Morozov, 

M. Roslova, O. Vakaliuk, S.-L. Drechsler, S. Johnston, D. 

V. Efremov, J. van den Brink, H. Rosner, M. Kumar, C. 

Hess, S. Wurmehl, A. U. B. Woher, B. Buechner, E. L. 

Green, J. Wosnitza, P. Vogt, A. Reifenberger, C. Enss, 

and R. Klingeler, arXiv::1301.5257. 
^® T. Shibauchi, private communication. 

A. V. Chubukov, Physica C 469, 640 (2009); S. Maiti, A. 

V. Chubukov Phys. Rev. B 82, 214515 (2010). 

A. J. Leggett, Prog, of Theor. Phys. 36 ,901 (1966). 
'^^ R. Cote, A. Griffin, Phys. Rev. B 48, 10404(1993); 

For the discussion on the plasmon mode in a dirty s-wave 

superconductor see B. Narozhny, A.I. Aleiner and B.L. Alt- 

shuler, Phys. Rev. B. 60, 7213 (1999); A. Kamenev, Field 

Theory of Non-Equilibrium systems, Cambridge University 

Press, 2011 and references therein. 

S.N. Artemenko and A.F. Volkov, Usp. Fiz. Nauk. 128, 3 
(1979) (Sov. Phys. Usp. 22, 295 (1979)); I.O. Kulik, O. 
Entin-Wohlman, and R. Orbach J. Low Temp. Phys. 43 
591 (1981); J.E. Mooij and G. Schon, Phys. Rev. Lett. 55, 
114 (1985); S.G. Sharapov, V.P. Gusynin and H. Beck, 
Euro. Phys. Jour. B 30 45, (2002); S.G. Sharapov and H. 
Beck, Phys. Rev. B 65, 134516 (2002). 
L. Benfatto, private communication. 

L. Benfatto, A. Toschi, and S. Caprara, Phys. Rev. B 69, 
184510 (2004). 
■'^ P.W. Anderson, Phys. Rev. B 112, 1900 (1958). 

Y. Ohashi and S. Takada, Journal, of Phys. Soc. Jap., 67, 
551 (1998) 

''^ L. Fanfarillo, L. Benfatto, S. Caprara, C. Castellani, and 
M. Grilh, Phys. Rev. B 79, 172508 (2009). 
R.V. Carlson and A.M. Goldman, Phys. Rev. Lett. 34, 11 
(1975). 



